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ABSTRACT 

Aims. The aim of the study is to constrain the physics of gamma-ray bursts (GRBs) by analysing the multi-wavelength 
afterglow data set of GRB 121024A that covers the full range from radio to X-rays. 

Methods. Using multi-epoch broad-band observations of the GRB 121024A afterglow, we measured the three charac¬ 
teristic break frequencies of the synchrotron spectrum. We used six epochs of combined XRT and GROND data to 
constrain the temporal slopes, the dust extinction, the X-ray absorption, and the spectral slope with high accuracy. 
Two more epochs of combined data from XRT, GROND, APEX, CARMA, and EVLA were used to set constraints on 
the break frequencies and therefore on the micro-physical and dynamical parameters. 

Results. The XRT and GROND light curves show a simultaneous and achromatic break at around 49 ks. As a result, 
the crossing of the synchrotron cooling break is no suitable explanation for the break in the light curve. The multi¬ 
wavelength data allow us to test two plausible scenarios explaining the break: a jet break, and the end of energy injection. 
The jet-break scenario requires a hard electron spectrum, a very low cooling break frequency, and a non-spreading jet. 
The energy injection avoids these problems, but requires e e > 1 (k — 2), spherical outflow, and cb < 10 -9 . 

Conclusions. In light of the extreme microphysical parameters required by the energy-injection model, we favour a jet- 
break scenario where < v sa to explain the observations. This scenario gives physically meaningful microphysical 
parameters, and it also naturally explains the reported detection of linear and circular polarisation. 

Key words. X-rays: bursts, gamma-ray burst: general, gamma-ray burst: individual - GRB 121024A, radiation mecha¬ 
nisms: non-thermal, stars: jets, methods: observational. 


1. Introduction 

Gamma-ray bursts (GRBs) are the most luminous phe^ 
nome na detected so far in the Universe ( |Klebesadel et al~ 
1973]) . They consist of pulses of gamma-rays emitted 


in a short time interval (milliseconds to hours) with an 
isotropic equivalent energy release of 10 5 ° — 10 54 erg. This 
is followed by fading multi-wavelength emission (fro m 


X-rays to radio), known as a GRB afterglow (Costa et al. 


1997} van Paradijs et al. 1997). In the standard afterglow 


model, the dominant process during the afterglow phase 
is synchrotron emission from shock-accelerated electrons 


Alexander von Humbolt Fellow 


in a collimated relativistic blast wave interacting with the 

k This unique 

I 


external medium 
featureless intrinsic spectrum 


(Meszaros & Rees 

makes' 


1997) 


B afterglows 

perfect events to study physical processes under extreme 
conditions (e.g. Fermi acceleration). 


The observed synchrotron spectrum is composed of four 
power-law segments joined at three main break frequencies. 
Each break frequency yields specific and correlated con¬ 
straints on the acceleration processes in the shock region 
and the dynamics and geometry of the relativistic outflow. 
Assuming an initially self-similar evol ution for the relativis- 
tic blast wave in the afterglow stage (Blandford & McKee 
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1976), the properties of the synchrotron spectrum can be 
expr essed in terms of constraints on the model parameters 
( Wijers et al.||1997 ; Sari et ah||1998 ). The dynamics of the 
outflow are dictated by the isotropic equivalent energy in 
the afterglow phase Ti SO , the circumburst medium density 
n = Ar~ k (where A is a scale factor, r is the radial 
distance from the source, and k is the power-law slope with 
k = 0 or k = 2 for inter-stellar medium (ISM) or stellar 
wind-like medium-density profiles, respectively) and, when 
the jet nature of the outflow becomes apparent, the jet 
half-opening angle #o (Rho ads||1999j ). The micro-physics of 
this synchrotron emission can be captured in a simplified 
manner using the post-shock energy fraction in accelerated 
electrons e e , the energy fraction in the magnetic field 6 b 
, and the power-law index p of the non-thermal electron 
population. 


A snapshot of the spectral energy distribution (SED) 
covering all break frequencies gives a number of ob¬ 
servational constraints equal to the number of model 
parameters. This scenario presents an ideal case in terms 
of model simplicity and data availability. In practice, 
the sample of bursts for which all spectral breaks could 
be simultaneously determined i s still small (about five, 
e .g. Panaitescu & Kuma r] 12002 ; Yost et ah] |2Q03| |Resmi| 


|010|. This 


et al. 2005 Chandra et al. [2008 Cenko et al. 
indicates the importance, in terms of sample statistics, 
of increasing the sample of bursts with simultaneous 
afterglow detections across the broadband spectrum. 
Common solutions when not all spectral breaks can be 
probed simultaneously are to extrapolate light curves in 
time and/or to fix one or more of the mo del parameters. 


Motivated by an apparent c orrelation (Frail et al. 2001 
Panaitescu & Kumar 2002), Ei SO has been equated to 
the total isotrop ic equivalent ene r gy release in gamm a 


rays E 1 ^ iso (e.g. Dai & Lu 1999 Pandey et al. 2003). 


The magnetic field energy densit y, e b has been fixed to a 


standard value (e.g. Frail et al. 2000 Cenko et al. 2010), 


or linked to the shock-accelerated electron en ergy, e e (e.g. 


[Medvedev 2006[ |van Eerten & Wijers 20091) . A standard 
valu e for the circumburst density has also been assumed 
(e.g. Chevalier et al.||2004 ). In these cases, the implications 
of the derived model parameters are conditional on the 
additional assumption(s). 

Here, we present the analysis of the simultaneous 
multi-wavelength observations of GRB 121024A. It was 
detected with the Swift satellite and had a redshift z 
= 2.30 measured with the X-shoo ter spectrograph at 
the Very Large Telescope (VLT) (jTanvir et al. 2012). It 
was followed up by different instruments in the radio to 
the X-ray regimes over several days. Linear and circular 
optical polarisation observat ions of the afterglow were 


taken (Wiersema et al. 2014). We report on the analysis 
of the broad-band SED of this afterglow, including X-ray, 
optical/NIR, sub-mm, and radio data. From these simul¬ 
taneous broad-band observations, we derive constraints on 
the micro-physical and dynamical parameters of the GRB 
afterglow. 


We provide a brief summary of the observations and 
relevant details of the data reduction in Sect. [2j We then 
describe in Sect. [3] a model-independent analysis of the 
data, starting with the description of the X-ray and op¬ 


tical/NIR light curves, followed by the description of the 
SED including effects of dust extinction and gas absorption 
at these wavelengths. In Sect. [4] we include radio and sub- 
mm data to study the broad-band SED in the framework 
of jet break and energy-injection scenarios. We derive all 
the micro-physical and dynamical parameters based on the 
standard afterglow model assumptions. In Sect. [5] we dis¬ 
cuss our results and compare the different viable scenarios. 
Finally, we conclude and summarise our results in Sect. [6] 


2. Observations and data reduction 


2.1. Swift 

On 2012 October 24 at To = 02:56:12 UT, the Swift Burst 
Alert Telescope (BAT, Barth elmy et al.|2005 ) triggered and 
located GRB 121024A flPagani et al. ]2012[ ). Swift slewed 


immediately to the bu rst, and the o bservations with the X 


Ray Telescope (XRT, Burrows et al. 2005) started 93 sec 
after the trigger. The observations were made in windowed 
timing (WT) mode during the first 242 s an d then w ere car - 
ried out in photon counting (PC) mode (P age et al.| 2012). 
The initial flux in the 0.2 — 10 keV band was 1.1 x 10 -9 
erg cm -2 s _1 . The Swift/XRT light curve and spectral 


data were o btained from the XRT repository (Evans et al. 
2007, 2009). The afterglow was located RA, Dec (J 2000) = 
04:41:53.28, -12:17:26.8 with an uncertainty of 078 (Pagani 
et al.||20l2| ) by the S wift/ UVOT, with a mag nitude in the 
b band of 18.4 db 0.2 (Holland & Pagani 2012). 


2.2. GROND 


The Gamma-Ray burst Optical Near-infrared Detector - 
GROND (Grei ner et al.|2 008) mounted at the Max-Planck- 
Gesellschaft (MPG) 2.2m telescope located at ESO La Silla 
observatory, Chile, was designed as a GRB follow-up in¬ 
strument. It provides simultaneous data in seven bands 
in a wavelength range from 400-2400 nm (g'r'i'z'JHK S ). 
GROND observations started 2.96 hours after the Swift 

and continued for the next 3.8 


2012 ) 


trigger (Knn st et al _ 

hours during the first night. The afterglow was detected 
in all seven bands at the position RA, Dec (J2000) = 
04:41:53.30, -12:17:26.5 with an uncertainty of 074 in each 
coordinate (Fig. [l]). After the observations during the first 
night, imaging of the field of GRB 121024A continued on 
the second, third, fourth, sixteenth, and seventeenth night 
after the burst. The optical/NIR data were reduced using 
standard IRAF tasks ( Tody|1993 Kriihler et al.|2008| ). The 
data were corre cted for Galactic foreground reddening E(B- 
V ) = 0.09 mag (Schlafly & Finkbeiner 2011), corresponding 
to an extinction of A v = 0.27 mag for R v = 3.1. The opti¬ 
cal magnitudes were calibrated against secondary stars in 
the GRB field (Table fl|. On 2 013 December 8 a Sl oan Dig¬ 
ital Sky Survey (SDSS) field (Aihara et al. 2011) at RA, 
Dec (J2000) = 04:59:42.0, -04:54:00 and the field of GRB 
121024A were consecutively observed in photometric con¬ 
ditions. The calibration of the secondary stars was made 
against the corrected zeropoints of the GRB field based 
on the SDSS field. The NIR magnitudes were cali b rate d 
again st the Two Micron All-Sky Survey (2MASS, Skrutskie 
et al][2006) catalogue stars in the field of the GRB. 
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Table 1: Secondary stars for photometric calibration. See Fig. [I] 


Star RA, Dec J(2000) £'(mag AB ) r'(mag AB ) i'(mag AB ) z'( mag AB ) J(mag Vega ) i7(mag Vega ) iF s (mag Vega ) 


I 

II 

III 

IV 

V 

VI 


04:41:49.55, 

04:41:52.36, 

04:41:55.40, 

04:41:47.91, 

04:41:46.29, 

04:41:52.83, 


-12:16:47.2 

-12:16:49.9 

-12:16:30.3 

-12:16:15.2 

-12:17:57.5 

-12:18:46.8 


19.96T0.05 

17.83i0.05 

20.73i0.06 

20.74i0.06 

17.76i0.05 

20.57i0.05 


18.75i0.05 

17.13i0.05 

19.45i0.05 

19.23i0.05 

17.36i0.05 

19.07i0.05 


18.21i0.06 

16.84i0.06 

18.91i0.06 

18.44i0.06 

17.21i0.06 

18.43i0.06 


17.92i0.06 

16.69i0.06 

18.68i0.06 

18.04i0.06 

17.11i0.06 

18.11i0.06 


16.73i0.06 

15.67i0.06 

17.45i0.07 

16.79i0.06 

16.19i0.06 

16.84i0.06 


15.96i0.07 

15.07i0.06 

16.66i0.07 

16.07i0.06 

15.81i0.06 

16.29i0.06 


15.88i0.08 

15.08i0.07 

15.96i0.06 

15.82i0.06 



Fig. 1: GROND r'-band finding chart. The secondary stars 
used for the calibration are labelled I-VI and are reported 
in Table [lj North is up and east to the left. 


2.3. APEX 


On 2012 October 24, we tri ggered an observation on the 
LABOCA bolometer camera (Siringo, G. et al. 2009) 1 Two 
observations at a frequency of 345 GHz with a bandwidth 
of 60 GHz were performed. The first started 19.8 ks af¬ 
ter the GRB, the second 98.7 ks after the GRB. On both 
days, the observ ations were taken in mapping mode and 
in on-off mode (Schuller et al. 112010). The data were re- 
duced using t he Bolometer Array analysis software (BoA, 
Sch uller| [ 2012 ). All the subscans (ten per scan) were used. 
A clipping of 2 a was used to remove any background ef¬ 
fects. The flux calibration was made using Jupiter for the 
focus, N2071IR as a secondary calibrator, and J0423-013 as 
a pointing source. There was no detection on either night, 
the upper limits are given in Table [2] 


2.4. Millimetre and radio observations 

In addition to the X-ray, GROND, and APEX data re¬ 
ported above, we also incorporated the following millime- 


1 Based on observations collected during Max-Planck Society 
time at the Atacama Pathfinder Experiment (APEX) under pro¬ 
posal m-090.f-0005-2012. 


Table 2: la Upper limits of the on-off measurements using 
the LABOCA instrument on APEX. 


Date 

On+off time [UTC] 

UL-Flux [mJy/beam] 

24-10-2012 

08:22 - 09:20 

3.6 

25-10-2012 

06:16 - 06:38 

10.4 

Date 

Mapping time [UTC] 

UL-Flux [mJy/beam] 

24-10-2012 

09:30-10:27, 10:39-11:00 

9.0 

25-10-2012 

08:52-09:29, 09:40-10:34, 
10:41-10:51 

12.0 


tre and radio observations reported in the literature in our 
SED analysis: 


The Combined Array for Research in Millimetre-Wave 
Astronomy (CARMA) started observations of the field of 
GRB 121024A ~ 120.9 ks after the BAT trigger at a mean 
frequency of ~ 85 GHz (3mm) ( jZauderer et al.| 2012). A 
mm counterpart was detected with a flux of 1.0 ± 0.3 mJy. 

The Very Large Array (VLA) started observations of 
the field of GRB 121024A ~ 109.0 ks after the trigger. A 
radio counterpart with flux of 0.1 0T0.03 mJy was detected 
at a frequency of 22 GHz (Laskar ^st al.|| 2012). 


3. Phenomenological data analysis 


We start with a model-independent analysis of the data. 
The observed flux is described by F ~ with a and 

/3 the temporal and spectral slope, respectively. First, we 
analyse the temporal evolution of the GRB 121024A af¬ 
terglow. Using the X-ray and optical/NIR light curves, we 
measured the temporal slope a and obtained information 
about particular features such as flares, breaks in the light 
curve, flattening, or any behaviour different from that ex¬ 


pected for a canonical afterglow light curve (LC) (Nousek 
Zhang et al. 2006). Then, we analysed the 
■ray to optical/NIR wavelengths at six differ¬ 
ent epochs. We obtained the spectral slope /3 and checked 
for spectral evolution. Given that absorption and dust ex¬ 
tinction only affect the data at X-ray and optical wave¬ 
lengths, we used this SED analysis to derive the host X- 
ray absorbing column density (V B ost ), which is commonly 
quoted as an equivalent neutral hydrogen column density, 
and the host visual dust extinction along the GRB line of 
sight (A^ ost ). 
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3.1. Afterglow light-curve fitting 


The temporal evolution of the X-ray afterglow of GRB 
121024A£]shows an initial steep decay with a temporal slope 
a = 3.6, followed by a small flare at ^ 300 s. For the present 
work, we only used the data after 10 4 s, which is the start 
time of our GROND observations (Figj2|. We tried to fit 
two models: First, a simple power law with host contribu¬ 
tion in the optical bands (plh) and slope decay a. Second, 
a smoothly broken power law with constant host contribu¬ 
tion ( brplh ) (EqJTj) (|Beuermann et al.||1999), with a pre and 
a pos t being the power-law slopes before and after the break, 
respectively, sm the smoothness, and t b the break time in 
the LC: 


( / _, \ -apresm / \ -a post sm'l 1 / srn 

Fu(t)=Cx |(£) +(£) ) +host. (1) 

The best fit to the X-ray light curve is a smoothly broken 
power law with a statistical significance x 2 /d.o.f. = 51/44 
(simple power-law: x 2 /d.o.f. = 87/47). The best-fit param¬ 
eters are an initial decay with <a pre = 0.84 ±0.09 and break 
time ££ rt = 32.5 ± 16.1 ks with sm = 5.0 db 2.6, followed by 
a steeper decay with <a pos t = 1.67 ± 0.23. 

The optical/NIR light curves (Table [3]) are well fitted 
by both a plh and a brplh model. A plh model gives 
X 2 /d.o.f. = 140/112 and a decay slope a = 1.07 ± 0.02, 
while a brplh model gives x 2 /d.o.f. = 107/109 and best¬ 
fitting parameters <a pre = 0.71 db 0.03, a post = 1.46 db 0.04, 
break time t£ pt = 31.4 ± 9.4 ks, and smoothness 2.7 ±1.1. 
Colour evolution in the optical bands is detected in the 
last epoch of our observations, which we associate with 
the increased contribution from the host galaxy. An F-test 
between the two model gives a null hypothesis probability 
of 1.86 x 10 -6 . Therefore, we conclude that the brplh 
profile describes the data in a better way. We used this 
profile for the following analysis. 

Both the X-ray and optical/NIR light curves are best 
fitted by a broken power law with similar break times. We 
therefore tried a combined fit to both the XRT and GROND 
light curves to test whether the same model can describe 
both data sets, which would thus better constrain the best- 
fit parameters. The best-fit model provides a good fit (x 2 = 
157 and d.o.f. = 141), with a pre-break temporal slope 
Q: p re = 0.86±0.05, post-break temporal slope <a post = 1.47± 
0.03, smoothness sm = 1.7 ± 0.3, and break time t\> = 
49.8 ±5.1 ks. 

3.2. Afterglow SED fitting 

We analysed six different spectral epochs using XRT and 
GROND data, spanning the time interval To ± 10 ks to 
Tq ± 240 ks, four before the break time in the light curve 
and two after it (Fig. p|. The spectral analysis includes the 
effect of the dust andmetal attenuation along the line of 
sight towards the source. For the last two SEDs, given the 
low signal-to-noise ratio at X-ray energies, we extracted the 
spectrum from the same time interval (40 ks - 240 ks), dur¬ 
ing which time there was no evidence of spectral evolution 
within the X-ray energy range. We then renormalised the 
spectra so that they corresponded to the measured X-ray 

2 http://www.swift.ac.uk/xrt_curves/536580 


flux of the afterglow at the mid-time of the two corre¬ 
sponding SEDs (i.e. tsED V = 88 ks and £sed VI = 107 ks). 



Fig. 3: Spectral energy distribution for the six SEDs high¬ 
lighted in Fig. [2| SEDs I - IV are from data before the 
observed break in the light curve. SEDs V and VI are from 
data taken after the break. The SEDs are scaled with an 
arbitrary factor for clarity in the plot. The values of /3 writ¬ 
ten above each line corresponds to the single power-law fit, 
where the slopes were left free to vary. The single power-law 
fit with a single tied slope has f3 = 0.86 ± 0.02. 


The SED analysis for all the six SEDs was performed 
simultaneously. The Galactic reddening was fixed to 
E(B — V ) = 0.09 mag,_corresponding_to an extinction of 
A& al = o.27 mag (Schlafly & Finkbeiner||2011) for a Milky 

Way (MW) reddening law. The Galactic absorbin g column 
~ ’ 20 , 1TY , —2 


density ± v ^ 


V^ al was fixed to 7.9 x 10 2U cm 2 (iKalberla 


et al.| |2005). The host magnitude was subtracted from the 
opticcd/NIR data, and the g' band was not included in the 
fit because of a damped Ly<a syst em along the line of sight 
towards the GRB (DLA, [FTiis et al. 2015). The values for 
the host extinction and absorption were tied between all 
the epochs, and the spectral slopes were left free to vary. A 
single power-law fit has a goodness of fit V7 d.o.f. = 28/46 
and all the spectral slopes values (see Fig.pJ) are consistent 
within 1 cr, confirming the lack of spectral evolution. A 
broken power-law fit either places the break outside the 
optical - X-ray frequency interval or fails to improve the 
fit when the break is forced to lie within this interval. In 
the latter case, the best-fit models have a goodness of fit 
x 2 /d.o.f = 32/40 


Given the lack of spectral evolution detected in our 
combined GROND/XRT light curve analysis out to 240 
ks, we fitted all six SEDs simultaneously with the same 
spectral model, with only the normalisation allowed to 
differ between epochs. The best-fit results are given by 
a single power law with a spectral slope /3 = 0.86 ± 0.02 
and goodness of fit x 2 /d.o.f. = 48/51. The best-fit 
host dust extinction given by jr Small Magellanic Cloud 
(SMC) reddening law ([Pd 1992) is A^ ost = 0.18 ± 0.04 
mag, and the host galaxy X-ray absorbing column is 
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Fig. 2: Light curve of the afterglow of GRB 121024A. Top: XRT light curve from the XRT repository. Bottom: GROND 
light curve in g'r'i'z' JHK S . The best fit for the combined light curve (optical/NIR and X-ray data) is a smoothly broken 
power law with host contribution, shown with dashed lines. The epochs used for the spectral analysis are highlighted 
with the vertical bars. The break time t\> = 49.8 i 5.1 ks. 


Table 3: Observed magnitudes of the GRB 121024A afterglow for the six highlighted epochs in Fig|2] The host contribution 
was subtracted. The magnitudes are not corrected for Galactic foreground extinction A^ al =0.27 mag. 


SED 

mid-time [s 

] ^(hiab) 

r'(m AB ) 

*'(m AB ) 

^(m AB ) 

J (my e g a ) 

H (myega) 

^s( m Vega) 

I 

11085 

20.75T0.08 

19.82i0.06 

19.53i0.06 

19.24i0.05 

18.68i0.10 

18.31i0.ll 

17.91i0.13 

II 

15497 

21.12i0.05 

20.20i0.05 

19.84i0.05 

19.63i0.04 

19.05i0.10 

18.67i0.10 

18.34i0.12 

III 

17006 

21.23i0.05 

20.32i0.05 

19.94i0.05 

19.68i0.04 

19.08i0.09 

18.68i0.10 

18.54i0.12 

IV 

21430 

21.48±0.27 

20.54i0.06 

20.20i0.06 

19.95i0.09 

19.43i0.10 

18.89i0.12 

18.81i0.15 

V 

88010 

22.89i0.24 

22.03i0.09 

21.74i0.10 

21.49i0.14 

21.01i0.26 

20.76i0.31 

20.41i0.31 

VI 

106998 

23.41i0.12 

22.54i0.08 

22.13i0.09 

21.85i0.13 

21.61i0.31 

21.18i0.33 

20.61i0.34 


N^ ost _ 0.30^0.29 x 10 22 cm -2 . The lack of spectral 
evolution together with the achromatic break observed in 
the light curve rule out the movement of the cooling break 
through the observed wavelength range. 

We extend the SED analysis in the following section 
with two additional epochs containing simultaneous obser¬ 
vations with XRT, GROND, APEX, CARMA, and EVLA. 


The first SED at t = 21.9 ks has GROND, XRT, and APEX 
data. The additional APEX upper limit requires a break 
between the APEX and NIR energies (see Fig. [4|. Then, 
we have a second SED at t = 109.0 ks, with two addi¬ 
tional measurements: CARMA and EVLA data points. The 
CARMA data point requires a break between the millime¬ 
tre and the NIR bands, and the EVLA data point implies 
a break between the radio and the CARMA wavelength. 
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Therefore at least two breaks in the broad-band spectrum 
of GRB 121024A are needed. These breaks are analysed in 
more detail in Sect. [4] in the context of the afterglow syn¬ 
chrotron spectrum model, where we use the constraints on 
/3, E(B — V) and N^ ost found in his section. 


4. Physical parameters of the standard 
afterglow model 

We now proceed with the derivation of the microphysical 
and dynamical parameters of the GRB afterglow, based on 
the standard afterglow model. In this model, the dominant 
emission is generally associated with synchrotron radiation 
from shock-accelerated electrons. These electrons are 
assumed to have a power-law energy distribution with 
slope p and minimum energy 7 m . The observed synchrotron 
spectrum is characterised by three main break frequencies 
(z/ c , Z'sa) and a peak flux. The synchrotron injection 
frequency is defined by 7 m . The cooling frequency v c is 
defined by the critical value y c , above which electrons ra¬ 
diate their energy on timescales shorter than the explosion 
timescale. The self-absorption frequency u sa marks the 
frequency below which the optical depth to synchrotron- 
self absorption is > 1. In this model, two main cooling 
regimes are defined by the relative position of the break 
frequencies: a fast cooling regime where z/ m > v c and most 
of the electron are cooling fast, and a slow cooling regime 
where < v c and most of the accelerated electrons are 


cooling slowly (Meszaros & Rees 1997 Granot & Sari 2002 ). 


The number of combinations of a and /? is limited when 
a specific dynamical model and the synchrotron spectrum 
are given. This gives rise to a unique set of relati ons be¬ 
tween a and {3_ known as M closure relations ” ([Meszaros fc 


Reesjl997 Sari et al.|1998 Zhang fc Meszaros. |2004p . These 

relations constrain the cooling regime, the circumburst en¬ 


vironment, the jet geometry, and the electron energy distri¬ 
bution p. We follow two main steps to analyse the afterglow 
data: 


1 . Spectral regime: The derivation of the p value and 
identification of the external density profile depend on 
the power-law segment of the synchrotron spectrum 
containin g the observing frequ ency. Using the closure 
relations (Racusin et al. 2009) together with the mea¬ 
sured parameters for a and 7 we find that the afterglow 
data can be described by two different spectral regimes 
(see Table [4]). In one case we have a spectral regime 
where v c < vk s (he. below the K s band), and in the 
other case we have one where v c > z/ xrt (i.e. above the 
XRT band). 


2. Microphysical and dynamical parameters: We include 
the APEX, CARMA, and EVLA data in our analysis 
and fit the data using a single, a double, or a triple 
broken power-law model depending on each individual 
case (in the double and triple broken power-law fits, we 
only consider sharp breaks because the data at millime¬ 
tre and radio frequencies are insufficient to constrain 
an additional free parameter i.e., smoothness). We use 
the standard formalism for a spherical blast wave prop¬ 
agating into an external cold medium during the slow 


cooling regime to de rive all the micro-physical and dy¬ 


namical parameters (Meszaros & Rees 1997 Granot & 
Sari| 2 QQ 2 ~ ), and subsequently check for consistency with 


the slow or fast cooling transition times. 


Both spectral regimes, that is, v c < z/k s z/ c > z/xrt, are 
explained and analysed in detail in the following subsection. 
The former corresponds to the jet in the light curve being 
associated with a jet break without energy injection, and 
the latter corresponds to the jet in the light curve associated 
with either the end of energy injection into the outflow or 
with a jet break with an ongoing energy injection during 
the whole evolution of the afterglow (i.e. the ongoing energy 
injection is still visible until the last observations). 


4.1. v c < vk s \ Jet break. 

Using the closure relations for a decelerating spherical blast 
wave, we find that the measured temporal slope before the 
break in the light curve is consistent with v c < vk s for both 
ISM and wind environments. This implies p— 1.73 ± 0.03, 
as (3 = p/2. The only plausible scenario consistent with 
the measured <a p0 st and /3 corresponds to a non-spreading 
uniform jet propagating into a wind environment. We 
therefore associate the achr omatic break observed in the 
light curve with a jet break ([Rhoads 1999 ; |Wiersema et al. 
2014). 


We proceed by including the post-break sub-mm and 
radio data in our analysis. The first broadband SED con¬ 
tains GROND, XRT, and APEX data. The best fit to this 
is a broken power law with both Galactic and host extinc¬ 
tion and absorption, with y 2 /d.o.f. = 3.6/5 (see Table |5j. 
The measured value of u c = 1.5 x 10 12 Hz is a lower limit 
because the APEX measurement is an upper limit. The 
second broadband SED contains XRT, GROND, CARMA, 
and EVLA detections. Two possible spectral sub-regimes in 
the slow cooling phase give a good fit to the data: The cool¬ 
ing regime where u sa < z/ m < z/ c , and the one where z/ m < 
z/ sa < z/ c - Because there are only a few data points at radio 
wavelengths, it is difficult to distinguish between these two 
cooling regimes. Therefore, we analysed both cases. 


Table 5: Results from SED fits for both a jet-break model 
with z/ sa < z/ m and z/ m < i/ sa , and for an energy-injection 
model. 


V 

Time 

Jet break 

Jet break 

Energy 

[Hz] 

[ks] 

Z'sa ^ Z7q 

Z'm^Z/ga 

injection 


26 

> 1.5 x 10 12 

> 1.5 x 10 12 

> 1.2 x 10 18 

109 

1.974 x 1012 

3.9 tli x 10 12 

> 1.2 x 10 18 

Z'm 

26 

- 

- 

< 1.4 x 10 14 

109 

1.371 x 1011 

< 2.2 x 10 10 

5.17J x 10 11 

Z'sa 

26 

- 

- 

- 

109 

8.3tR x 10 10 

7.4 ±o .7 x 10 10 

7.47.6 x 1010 


When 1 < p < 2 (i.e. a hard electron spectrum), there 
is more energy-per-decade in high-energy electrons. This 
distribution has important implications for the analysis of 
the physics in the shock region, specifically requiring an 
additional high-energy cut-off in the electron population. 
We based our analysis on the assumption of a proportion¬ 
ality between 7 m and 7, where 7 is the Lorentz factor of 
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Table 4: Closure relations*. /3 = 0.86±0.02 was use in the analysis. When determining the energy-injection 
parameter g, we use the measured a pre = 0.86 ± 0.05 and <a pos t = 1.47 ± 0.0 3. The equatio ns used for q 
are for the case when p > 2, for both spectral regimes. More details in Racusin et ah||2009 


Instantaneous injection 

a(/3) 


Energy injection 
q(P,a) 


Section* 


Spherical outflow 


Vm < V Q 


ISM 

Wind 


3/3+5 

8 

/3+3 

4 


0.95 ± 0.01 
0.96 ± 0.01 


2(l+q—ft) 
ft +1 

2(l+a—ft) 
ft+l 


1.07 ±0.05 
1.07 ±0.05 


4.1 


Z+n < Z/ < V c 


ISM 

Wind 


3ft 

2 

3/3+1 

2 


1.29 ±0.03 2(1 +^~ /3) 


= 1.79 ±0.03 


2(a-ft) 

ft+1 


= 0.69 ± 0.04 

= -0.01 ±0.05 


4.2 


Uniform non-spreading jet 


^c/m ± ^ 

ISM 

Wind 

3/3+11 

8 

ft+5 

4 

= 1.70 + 0.01 

= 1.47 + 0.01 

2(l+2a—2ft) 

3+2/3 - 

2(1+0!—ft) 

2+ft 

0.95 ±0.05 

1.13 ±0.04 

4.1 


t'm < V < P c b 

ISM 

Wind 

6/3+3 

4 

3/3+2 

2 

= 1.64 ±0.01 

= 1.84 ±0.01 

2(l+2o—2ft) 

5+2/3 “ 

2(a-/3) _ 

2+ft 

0.65 ± 0.03 

0.43 ± 0.04 

4.3 


Uniform spreading jet 

^c/m ± ^ 

ISM/Wind 

ft+3 

2 

= 1.93 ±0.01 

2+3o-4ft 

2(/3+l) 

0.79 ±0.05 

4.3 


Vm < V < V c b 

ISM/Wind 

2/3+1 

= 2.72 ±0.04 

l+3o—4ft 

2(/3+2) 

0.35 ±0.03 

4.3 



When v c < v < Vm = 0.5 and it does not depend on p or a. We did not include this scenario as it is not 
compatible with our data at any time. 

Details on the results and implications of the closure relations are discussed in the outlined section. 
a For is c >is,p = 2 p. When p = 0.86 ± 0.02 we have p = 1.73 ± 0.03 (1 < p < 2). 
b For z+i < v < v c . When p = 0.86 ± 0.02 we have p = 2.73 ± 0.03 (p > 2). 


the shocked fluid. This implies that is proportional 
to local temperature, which is physically plausible since 
the non-thermal population is presumably accelerated 
out of a Maxwellian population. The upper cut-off in 
the electron distribution can be assumed to lie beyond 
the X-ray band and does not need to be accounted for 
explicitly. e e can no longer be interpreted as the fraction of 
energy in accelerated electrons. Instead, it becomes a scale 
factor between 7 a nd + m , according to 7 m = K x 7 with 




1.07 x 10 5 s when z/ m < z/ s; 


x m p /m e (Meszaros fc Rees| 


1997 


K 

the formalism used T5}7 liranot &; Sari (200! 


I, 


We followed 
who derived 

the flux equation usi ng a full fluid profile for the blast 
wave (Blandford fc McKee|[ l976) and took the line-of-sight 
effect and the cooling times oFThe 
account. 


re individual electrons into 


Based on the values for the break frequencies presented 
in Table [ 5 ] for both spectral regimes, z/ sa < z/ m (see Fig. 4a) 
and z+n < z/ sa (see Fig. 4b), we derived the microphysical 
and dynamical parameters. The results are reported in 
Table H and are used to calculate the transition times 
between the spectral regimes. First, the transition from 
fast to slow cooling. This corresponds to t Uc=Vrn ~ 2.8 x 10 4 
s and t Uc=Vrn ~ 2.6 x 10 3 s for z/ m > z/ sa and v m < z/ sa , 
respectively. In both cases, it is before the time of the 
analysed SED at t — 109 s, confirming the slow cooling 
assumption. Second, the transition from optically thin 
to optically thick, that is, when z/ m goes below z/ sa . This 
occurs at ^ 1.8 x 10 5 s when z/ m > z/ sa and at 


4.2. u c > zz X rt : Energy injection. 

The closure relations (Table. [4| and the possible spectral 
break positions resulting from fitting synchrotron spectra 
to the SED allow for an alternative scenario, where 
v c > vxrt and p > 2. In this case, the break between 
the mm and NIR wavelength corresponds to z/ m and the 
break in the light curve is associated with the end of the 
ongoing energy-injection phase. Smooth energy injection 
into the ejecta can result from slower shells with a range of 
velocities that catch up with each other or from a long-term 
engine luminosity. In the latter case, the energy-injection 
parameter q is defined by L = Lo(t/tb) -q . Using the flux 
and frequency eq uatio ns for radial flow from |van Eerten| 


& Wijers (2009) and Leventis et al. |2012 , we derived 


the closure relation for a general density profile with an 
arbitrary k during the deceleration stage following energy 
injection (fc = for < v < z/ c ). The best-fit 

results for a post and P then imply k = 1.05 ± 0.23. 

During the energy-injection phase, a forward-reverse 
shock system is set up in the flow. Using the flux equa¬ 
tion describing the energy-injection phase from van Eerten 


(2014), we have a given relation between <a, /?, fc and q. If 
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the emission is dominated by that from the forward shock 
(FS): 


FS : q = 


8 — 2a (—4 ± k) ± 2/3(—4 ± k) — Ak 
3fc-8 + /3(Jfe-4) 5 


( 2 ) 


and the following relation if the reverse shock (RS) emis¬ 
sion is dominant: 


RS 


Q = 


4 ± 8a — 2(1 Hff ' ol ± /3)/c 
(3 + /3)k — 10 ' 


(3) 


The values for <a pre and /3 derived in Sect. 3.1 imply 


q = 0.52±0.07 for dominant FS emission and q = 0.88±0.09 
for dominant RS emission. These q values (as well as the 
pre-break temporal slopes) are consistent with those deter¬ 
mined for Swift samples (e.g. |Racusin et ah]|2QQ9 Evans 
Margutti et al.|2013|. If we fix k = 2, we obtain 


et al. 2009; 
q — 0 for b 


oth RS a nd FiS, consistent with predictions for 


q ~ 

a magnetar model (Dai & Lu 


2004). For an ISM density profile, q 


1998 


Zhang & Meszaros. 
= 0.69 ± 0.04 for FS 
emission and q = 1.09 ± 0.03 for RS emission. RS emission 
can therefore not be dominant because q > 1 implies that 
the energy injection decays too rapidly to sustain a plateau. 


approximation leading to the closure relations used here 
is no longer valid. On balance, the bulk of the energy wil l 
remain confined to 0o (see discussion in van Eerten 2014), 
and then the non-spreading jet approximation is favoured. 
The non-spreading case gives q = 0.65 =b 0.03 for an ISM 
density profile and q = 0.43 ±0.03. Therefore, only the ISM 
density profile is consistent with the pre-break q value, and 
the prolonged energy injection would only be possible if 
the observed jet break is due to geometrical effects alone. 


We analysed our last two SEDs including radio, submm, 
NIR, optical, and X-ray data. We used the flux and break 
frequen cy equ ations for energy injection presented in |van| 
Eerten (2014) together with equations for u sa (van Eerten 
in prep.) to obtain the model parameter values presented 
in Table [6] As in the previous case, where energy injection 
was only operating before the break in the light curve, some 
unphysical values for the parameters are found. The main 
problems are e e > 7.6, when it should not be greater than 
unity, and the value for the density ~ 10 7 cm -3 , instead of 
being of order unity as expected. 


5. Discussion 


After the energy-injection phase, only a decelerat¬ 
ing forward shock remains, and a standard afterglow 
emission model can be applied. We therefore proceeded 
with the analysis of the final SED at t = 109s, which 
contains EVLA, CARMA, GROND, and XRT data. 
The best-fit profile is a sharp double broken power law 
with x 2 /d.o.f. = 8.50/8 (Fig. [4c]). The critical values 
reported in Table [5] were used to derive the micro-physical 
parameters after the energy-injection phase (see Table 
[6|. In this scenario, v c cannot be measured and we can 
only place a lower limit. The k value lies just between 
the expected values for ISM and wind environments, and 
therefore we determi ned the values for bo th wind and ISM 
envir onments using Granot & Sari ( 2002|) and for k = 1.05 
using ” " . 


van. 


Eerten & Wij ers (2009) and |Leventis et ah] ( |2012| ). 


4.3. v c > z/ xrt : Energy injection and jet break. 


Now, we analyse the afterglow parameters assuming 
prolonged energy injection at all times, both before and 
after the break in the light curve, and the break in the 
ligh t curve is associated with a jet break. As shown in Sect. 
|4.2[ RS emission is not dominant before the break, and 
therefore we assumed only a dominant FS emission. To 
avoid too many free parameters, we restricted the study 
to ISM and wind density profiles. The analysis for these 
two medium profiles for the ongoing energy-injection phase 
befo re the break in the light curve is presented in Sect. 
(4.2). Here, we used a pos t to derive the q values after the 
break and, assuming q does not evolve, we compared these 
post-break values with the pre-break values to check if it 
is possible to have ongoing energy injection together with 
a jet break. We analysed two cases for the post-jet-break 
evolution: a sideways spreading jet and a non-spreading jet. 
For the former case q = 0.35 ± 0.03, inconsistent with the q 
value before the jet break. In addition to this inconsistency, 
if the energy is continuously injected within 0 O? while the 
front of the jet begins to spread, the homogeneous shell 


In the previous sections we have presented a detailed anal¬ 
ysis of the afterglow observations and derived values for the 
microphysical and dynamical parameters. Here we compare 
the derived values in the different scenarios. We discuss the 
positive and negative aspects of each model in the frame¬ 
work of the standard afterglow model. 


5.1. Jet break without energy injection 


This scenario requires three main features: First, the 
cooling break must be at around a few times 10 12 Hz 
at > To ± 109 ks. Although a value for v c th is low has 
been seen before (i.e. GRB 060418, Cenko et al. 2010), 
in more than 95% of a combined GROND-XRT sam ple, 
v c was dete cted above the optical frequencies (Greiner 
et al. 2011). Second, the closure relations require that 


the jet does not spread out sideways following the break 
time. The jet has to remain in this non-spreading state at 
least until about one day after the jet break because no 
spectral evolution is detected so far in the observations. 
This behaviour is at odds with findings from theoretical 


iGranot & Piran 2012) and numerical (van Eerten et al. 

>0T0l 


|De Colic et al |2012 |van Eerten &; J\4 acFady< 


en 


2012) studies of afterglow jets. Third, a very hard electron 
spectrum (p < 2) with p = 1.73 is required, as are 
additional assumptions about the minimal Lorentz factor. 
Although this is significantly lower than the value_of 2.3 
expected from Fermi accele ration theory (e.g. Kirk et al. 


2000 

of vaJ 


Achterberg et al.|2 001), it is within the average range 
ues 1.5 — 3.0 observed in previous GRB afterglow 


studies (Curran et al. 2010). 


Two diff erent spectral sub-regimes were presented 
in Sect. |4.1[ either with z/ sa < i/ m or with z/ m < z/ sa . 
Here the main assumption is 7 m ^ 7. For both spectral 
sub-regimes, the derived values for cb (Table [6| are in 
the same range as previous measurements reported in 
the literature, and the values of 0q are consistent with 
a collimated outflow (0.1 — 0.3 rad). The values for the 
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and j4.3 e P = e^x (| p — 2| )/(p — 1) 


Granot et al. (2005). n = Ar 2 with 


Table 6: 7 m ,e e , 6 b, E- lso , n and 0q for the models described in Sects. |4.l[ |4.2 , 
and -E'i SO ,52 = Fiso/10 52 . The half-opening ang le is der ived using Eq.(4) from 

A = M/47rv w = 5 x 10 n A* g cm -1 (Chevalier & LJ2000). For k = 2 we report the density in terms of A*. For k = 0 and 
k = 1.05 we report the number density no in units of cm -3 . In the special case of k = 1.05 we use a reference distance 
of r = 10 17 cm. 



^B 

A * , n 0 

0o [rad] 

■E'iso,52 [erg] 


Jet break, p = 1.73 db 0.03, v c < z/ Ks 


^sa ^ 

102.7t5f 2 6 

2.09ij'o5 • 10-2 

2 - 117 : 9 ? • 10- 2 

1 41 + 4-01 

Q 02 + 0 - 07 
u * oz -0.02 

q r+ 0.07 
u * iO - 0.03 

987% 

Z^m ^ ^sa 

< 11.2 

< 9.31 • KT 4 

< 7.87 • 10“ 2 

> 0.78 

> 0.13 

> 2.94 

< 74% 


Energy injection until £5 

in the light curve, p = 2.73 ± 0.03, v c > z/ xrt 


II 

to 

> 2.01 • 10 3 

> 1.10 

< 6.64 • 10" 10 

> 1.23 • 10 3 

>0.85 

> 2.36 

< 78% 

k = 1.05 

> 1.4 • 10 3 

> 0.76 

< 2.1 • 10" 9 

> 4.3 • 10 5 

>0.8 

> 3.4 

< 71% 

O 

II 

> 1.11 • 10 3 

> 0.75 

< 2.25 • 10" 9 

> 1.21 • 10 7 

>0.77 

> 3.67 

< 69% 


Energy injection scenario with jet break, 

p = 2.73 d= 0.03, v c > zz xrt 


O 

II 

> 1.6 • 10 4 

> 7.6 

< 3.9 • 10 -8 

> 1.29 • 10 7 

> 1.21 • 10" 2 

> 0.16 

< 98% 


circumburst density therefore agree with the collapsar 
model and a Wolf-Rayet star as possible progenitor, with 
mass-loss rates of ~ 1.4 x 10 -5 M^yr -1 when z/ sa < 


and > 7.8 x 10 b M 0 yr 1 when v ru < n sa , for a wind 


velocity v = 1000 km s 1 (Chevalier & Li|1999[[2000| ). The 
efficiency requirements are extremely high, For z/ sa < i/ m 
the measured Ei SO implies an efficiency of p ~ 98% and for 
v m < ^sa the efficiency is r\ < 74%. These two efficiency 
values are much higher than expected in the standard 
fireball sh ock model, for which an efficiency of 77 < 10% is 
predicted (iKobayashi et al. 1119971 iDaigne & Mochkovitch 

- nariPT" .. . ' 1 M . 


1998 Kumar||1999[ Cranot et al.||2006[ t^enko et al.||2011| )T 


As a final verification of th is model, we applied a condi¬ 
tion from Granot & Sari ( 2002 ) to the evolution of the after¬ 
glow spectrum in a wind environment for a given set of mi¬ 
crophysical and environmental parameters. This states that 
if E is ^J 2 e^ 7 > 100 , the afterglow spectrum evolves 

from fast to slow cooling, where in the slow cooling phase, 
initially z/ sa < z/ m < z/ c , but eventually u m < z/ sa < v c . If 
A * e"" e “ 1 E^ 52 7 < 100 , the afterglow spectrum only goes 
through one spectral regime in the slow cooling phase where 
v m < * 4 a < v c . In this latter scenario we are never in the 
regime where z/ sa < during the slow cooling phase. In 
our jet-break model where the spectral regime is z/ sa < z/ m , 
we therefore require that the derived micro-physical and 
dynamical parameters give A^e^ 1 E is °> 100. How¬ 
ever, we find that our best-fit values presented in Table 
m for V a a < v m give 52, inconsistent with the condition 
stated above, and therefore this regime can be ruled out. 


3 Efficiency 
the outflow 


V — EisOj-y /(E 

iso ,7 


of the conversion of the kinetic energy in 
to gamma-rays during the prompt emission 


+ Ei S o). Fi so ,7 i s the isotropic energy 
released in the prompt gamma- rav emission. In this case 

171. _ O J+2.6 

-f-^iso,7 — 2.2 


x 10 52 erg (|Butler & Kocevski 2007) 
(http: //butler .lab.asu.edu/Swift / index. Jit ml ). it is calculated 
using Fi so ,7 — 47rdLF/(l + z), where F is the fluence in the 
gamma-ray band. BAT: from 15 — 150 keV in the observer-frame. 
Fi so : energy range 1 — 10 4 keV in the rest frame. 


The favoured regime is then a slow cooling phase where 
v m < where our best-fit parameters give the value 257. 


5.2. Energy injection 


According to the shape of the spectrum and the closure 
relations, it is also possible to have v c > z/x rt? implying an 
energy-injection model. In the energy-injection scenario, 
both with and without a jet break, our best-fit values for 
q are consistent with q ~ 0.5, corresponding to smooth 
energy injection, _which has been observed in several other 

electron 
have 


Zhang et al. 2006), and the hard 


cases (e.g. 

spectrum is not required anymore since we now 
p = 2/3 + 1 = 2.73. However, further problems with the 
other afterglow parameters are found. 


The energy-injection scenario together with a jet break 
in an ISM external medium and without a jet break in 
a wind-like (k = 2) external medium can be ruled out 
because e e > 1, therefore, this scenario is physically not 
meaningful. There are two other scenarios to be analysed 
then: the energy-inject ion scenario without a jet break in 
an ISM medium and with a general density profile with 
slow k = 1.05. The analysis below focuses on these two 
cases. 


In relation to the micro-physical parameters, the 
derived values for 6 b differ from previous observations, 
but agree with theoretical predictions. In the former 
case 6 b is more than four orders of magnitude (< 10 -9 ) 


studies (e.g. Panaitescu & Kumar 

2002; Yost et al. 

2003; 

Panaitescu 2005 Cenko et al. 2010 

). In the latter case, the 


with expected values from shock compression of the seed 
m agnetic field (B 0 ^ pG ) in the _surrounding medium 


(Kumar & Barniol Duran 2009 Santana et al. 2014), and 


no further amplification or additional magnetic field would 
be required in the shock region. On the other hand, the 
derived value for e e is < 0.75 consistent with theory as 
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e e rsj 0.2 (Sa ntana et al.| |2014). 





(c) Energy injection 

Fig. 4: Broadband SEDs of the afterglow of GRB 121024A 
from the radio to the X-ray regime for the three models de¬ 
scribed in Sect. [U Blue line: SED at t = 21.9 ks. Red line: 
SED at t = 109.0 ks. The dashed lines represent the ab¬ 
sorbed model, the solid lines the unabsorbed model. The 
grey-shaded regions corresponds to the la limits of the 
model. 


e e < 1 but higher than the average of observations where 


Furthermore, in the case where cb <C e e , as implied by 
our analysis, we would expect there to be a contribution 
to the cooling of el ectrons fro m i nverse C ompton (IC) 
scattering processes (Panaitescu & Kumar 2000 Sari & 
Esin 2001). The IC emission will mainly affect the cooling 
frequency in the slow cooling regime. The final value of v c 
is expected to be lower than the value with synchrotron 
cooling alone. The IC contribution to the total observed 
afterglow emission can be included using the Compton 
parameter defined as 7 = rj c e e /e b, where rj c = { r y c / lm) 2 ~ v 
for t he slo w cooling regime (for more details see |Sari & 
Esin 2001). With this parameter the cooling frequency 
will be lowered by a factor of (1 + Y) 2 . A constraint C 
can be derived to test whether the IC contribution is 
important during the evolution of the observed emission 
from the after glow, or if it can be neglected (see Eq. (4.9) 
of Sari & Esin 2001). This constraint depends only on 
the observationalquantities (break frequencies and peak 
flux) of the afterglow and is independent of the theoretical 
afterglow parameters. It can be expressed in terms of Y 
as C = Y/(l 4- Y) 2 and has a consistent solution for Y 
only if C < 1/4. Using our measurements, we find C ~ 10 6 
for ISM. This result indicates that the IC component is 
not a relevant contribution for this afterglow, contrary 
to the theoretical expectation when e^/e e « 1 and the 
energy-injection scenario is not favoured. 


The lower limits derived for the density are 2 (k = 1.1) 
and 4 (k = 0) orders of magnitude higher than previous 
density measurements for bursts with similar isotropic en- 
ergies ( E iso rsj 10 52 erg) between 10 2 and 10 3 for constant 
dens i ty cir cnmbnrst profiles (n) (e.g. see Fig. 11 in Cenko 


et al. 2011). The values for 9q for the case of no jet break 
indicate a spherical outflow, opposite to the collimated out¬ 
flow usually expected and assumed in the standard after¬ 
glow theory. Finally, the energy lower limit is ^ 3 x 10 52 
erg, implying an efficiency of r] < 70%. 


5.3. Origin of the light-curve break 

From the available data, it is difficult to make a clear case 
for a preferred model for this GRB afterglow. Each of the 
studied models has specific problems that are difficult to 
explain with a simple afterglow model and would probably 
be better understood with a more complex and detailed 
mod el of the afterglow emission, especially at early time s 
(e.g. Waxman & Meszaros 2003 Morsony et al. 20071 ). 
However, we are able to rule out some of the possible 
models. For instance, the jet-break model where the spec¬ 
tral regime is u sa < z/ m is ruled out because the spectral 
evolution will n ever cross that regime in the slow cooling 
phase (see Sect. 4.1). In a similar way, we can rule out the 


energy-injection model with a wind density profile k = 2, 
and the energy-injection model with a jet break for k = 0 
because e e has to be larger than one, which is physically 
not meaningful. The energy-inject ion model without jet 
break for k = 1.1 and ISM density profiles cannot be 
ruled out. However, in this model the extremely high 
density requirements are far from theoretically expected 
values and previous measurements. Moreover, the resulting 
spherical outflow geometry, implied by the derived value 
for the half-opening angle, would require a very energetic 
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explosion. 


For more than 40% of the X-ray afterglows an initial 
plateau last i ng for about 10 3 — IQ 4 s is observed ( Lazzati fc 
Perna||2007 Margutti et al. 2013) and has been associated 


with a co ntinuous energy inject i on during the afte rglow 


evolution ( |Nousek et al. |2006| van Eerten 
’ th< ’ ~ 


2014). We 


compared the X-ray luminosity (0.3 — 30 keV) (Lf) and the 
break time in th e rest frame (tf F ) to the relation observed 
in other GRBs (Dainotti et al 


1[2008). Figure |5] shows the 
Observed in a sample of 62 


relation between L~f and observed in a sample 
long GRBs studied by|Margutti et al. (2013|). We included 
GRB 121024A for both an energy-inject ion phase that ends 
at the time of the break in the light curve (red star) and for 
an ongoing injection phase until the end of the observations 
at t = 240 ks, both taking the end-time luminosity directly 
(green star) and correcting for the change in the slope 
introduced at the break (grey star). Energy injection 
follows the correlation very well, supporting the scenario 
with energy injection up to the break in the light curve. 
Continued energy injection is disfavoured in view of the 
correlation. We note that assuming energy injection to ex¬ 
tend beyond the final data point at 240 ks will only shift the 
grey and green stars even farther away from the correlation. 



Fig. 5: 0.3 — 30 keV luminosity at th e end of the energy- 
injection p hase dDainotti et al. 2008). The black dots are 


taken from Margutti et al. (|2013). The stars correspond to 
the GRB 121024A afterglow: the red star to when the end of 
the plateau phase lies at 49.8 ks, the green star corresponds 
to the ongoing energy-inject ion phase before and after the 
break in the LC at t = 240 ks, and the grey star corresponds 
to the luminosity corrected for the jet-break effect. The 
dashed line in the middle corresponds to the best fit, and 
the shaded region is the la error of the fit. 


We consider the jet-break model where z/ m < z/ sa to be 
the preferred scenario. In this model, all the micro-physical 
and dynamical parameters are within the range of previ¬ 
ous measurements and within the expected values from the 
standard afterglow model. The low values for the energy 
are just lower limits and therefore are not a strong ar¬ 
gument against this model. The main problem is related 
to the hard electron spectrum that requires additional as¬ 


sumptions on the acceleration process of the electrons in 
the shock region. However, this is certainly not the first 
GRB for which such a shallow electron spectrum was de¬ 
rived, and viable ways to handle this scenario have been 
put forward, two of which we investigated and found to 
give reasonable and physically meaningful results. The de¬ 
rived hard electron spectrum need not be a reason to reject 
a model, and more likely reflects our poor understanding 
of acceleration processes under extreme conditio ns. Fi nally, 
the li near polarisation observations reported by |Wie rsema 


et al. (2014) would agree with a jet-break model where the 
linear polarisation would be a direct result from the jet 
break. However, there are still no studies reported in the 
literature analysing whether it would be possible to obtain 
this type of polarisation from an energy-inject ion model. 


6. Summary and conclusions 

We analysed the afterglow of GRB 121024A and showed 
that the multi-wavelength data enable the afterglow 
spectra and temporal parameters to be measured to a high 
degree of accuracy, which sets strong constraints on the 
micro-physics in the shock region and on the dynamics of 
the jet. The combined GROND and XRT data allowed us 
to determine the spectral slope /3 in this energy regime 
with high accuracy, and therefore we were able to measure 
the electron index p. We modelled our complete set of 
observations using two different physical interpretations: 
a jet-break model and an energy-inject ion model. The 
energy-injection model requires 77 < 77%, 71%, and 69% 
for k = 2 , 1 . 1 , and 0 , respectively, and does not contradict 
Fermi acceleration predictions for the electron index p. 
However, it does face some problems with the derived 
microphysical parameters in the case of a wind density 
profile, and the density values are extremely high in all 
three of the density profiles studied. 

The jet-break model requires a hard electron spec¬ 
trum. Assuming oc 7 , the derived microphysical and 
dynamical parameters are all consistent with previous 
measurements and with expected values from theoretical 
analysis. There is a problem with the efficiency require¬ 
ments, which for z/ sa < u m can be as high as 77 ~ 98%, and 
for < z/ sa the efficiency has an upper limit of 77 < 71%. 

The results presented here on GRB 121024A show that 
broadband afterglow data from the X-ray to the radio al¬ 
low for a detailed analysis of the characteristic properties 
of the GRB afterglow synchrotron emission spectrum. As 
studies of other GRBs have also shown, such datasets are 
invaluable for determining the range of microphysical and 
dynamical parameters within GRB shock-fronts with bet¬ 
ter statistics and avoiding adding additional assumptions 
to the analysis. Through our extensive data coverage of 
GRB 121024A, we have been able to constrain the position 
of all synchrotron breaks, which in turn has allowed us to 
measure, or set constraints on, all the micro-physical and 
dynamical parameters of a GRB afterglow. This informa¬ 
tion is crucial to further study the GRB afterglow emission 
process, and we are currently working on a larger sample of 
GRBs with sensitive and broadband afterglow data (Varela 
et al. in preparation). Future continual coverage of the GRB 
afterglows with sensitive telescopes over a wide wavelength 
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range and at multiple epochs will enable us to place strong 
constraints on the micro-physical parameters for a larger 
sample of GRBs, and allow us for instance to investigate 
the evolution of these parameters. 
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Appendix A: Alternative analysis to the 
jet-break model 

One of the interpretations for the break in the light curve 
is a jet break. As mentioned in Sect. |4.1[ a hard electron 
spectrum is required in this scenario. For completeness, we 
present he re a second approach to analyse a hard electron 
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'ect of 


spectrum QBhattacharya| [2001| |Dai fe Cheng 
et al.|2013 ). Here, instead of assuming y m oc 7 , the 
an upper cut-off ym = [3e/0<TxB ] 1 / 2 in the energy range of 


IS 

eff< 


the accelerated particle population can be included in the 
minimal Loren t z facto r such that y m oc 
(Dai & Cheng 2001). This upper cut-off 


(77m 


(p-2)\(l/p-l) 
M ) 

follows from 

equating acceleration and synchrotron cooling timescales. 
The advantages of this approach are that the extra cut-off 
is modelled explicitly and that e e can still be interpreted as 
the fraction of energy in the accelerated electrons. However, 
it implies that the behaviour of the electron population 
at low energies is dictated by the few electrons that were 
accelerated more efficiently, which is not supported by any 
clear physical mechanism. 


We calculated all the parameters for both regimes, when 
z^a < and when z/ m < and report them in Table 
The transition from the fast to slow cooling regime 
occurs at t Uc=Vrn ~ 3.6 x 10 4 s for z/ sa < v m and t Vc = Vrn ~ 
6.6 x 10 2 s for z/ sa > consistent with the slow cooling 
assumption. The transition from optically thin to optically 
thick occurs at t Usa=l/ ~ 9.8 x 10 4 s for z/ sa < z/ m and 
at t Vsa=Vrn ~ 1.2 x 1(? s for z/ sa > z/ m . In this case, the 
sub-regime where z/ sa < z/ m is ruled out because here the 
time that z/ m would cross z/ sa is ^ 12 ks, which is before 
the epoch of the SED (109 ks) used in the analysis. In the 
second sub-regime, where z/ m < z/ sa , 0 q is consistent with a 


collimated outflow and A* is in the range of expected values 
for a wind environment and corresponds to a mass-loss rate 


-1 


of >2.6x10 6 M 0 yr 1 for a wind velocity v = 1000 km s 
(Chevalier & Lil 1999 2000), consistent with a Wolf-Rayet 


star as a possible progenitor. The efficiency rj < 98%, even 
though it is just an upper limit, is extre mely high. Our 
results agree with the results in Sect. |4.1| and support the 
jet-break scenario when v m < v sa as the preferred scenario. 


Table A.l: y m ,e e , e B, ^iso> n and 60 for a jet-break model. 
p = 1.73 ± 0.03. z/ c < vk s 


C e 

cb 


^sa ^ 

Z^m ^ ^sa 

0-80™ 

< 0.11 

7-46te j? • 10- 3 
< 0.18 

2-07™ 

> 0.26 

60 [rad] 

E iso , 5 2 [erg] 

f ] 

^sa ^ 

Z^m ^ ^sa 

3.75™ • 10- 2 
> 1.12- 10“ 2 

1.25^0.75 • 10 -2 
> 0.19 

99t|% 

< 98% 


References 

Achterberg, A., Gallant, Y. A., Kirk, J. G., & Guthmann, A. W. 2001, 
MNRAS, 328, 393 

Aihara, H., Allende Prieto, C., An, D., et al. 2011, ApJS, 193, 29 
Barthelmy, S. D., Barbier, L. M., Cummings, J. R., et al. 2005, 
Space Sci. Rev., 120, 143 

Beuermann, K., Hessman, F. V., Reinsch, K., et al. 1999, A&A, 352, 
L26 

Bhattacharya, D. 2001, Bulletin of the Astronomical Society of India, 
29, 107 

Blandford, R. D. & McKee, C. F. 1976, Physics of Fluids, 19, 1130 
Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2005, Space Sci. Rev., 
120, 165 

Butler, N. R. & Kocevski, D. 2007, ApJ, 663, 407 
Cenko, S. B., Frail, D. A., Harrison, F. A., et al. 2011, ApJ, 732, 29 

Cenko, S. B., Frail, D. A., Harrison, F. A., et al. 2010, ApJ, 711, 641 

Chandra, P., Cenko, S. B., Frail, D. A., et al. 2008, ApJ, 683, 924 
Chevalier, R. A. & Li, Z.-Y. 1999, ApJ, 520, L29 

Chevalier, R. A. & Li, Z.-Y. 2000, ApJ, 536, 195 

Chevalier, R. A., Li, Z.-Y., &; Fransson, C. 2004, ApJ, 606, 369 
Costa, E., Frontera, F., Heise, J., et al. 1997, Nature, 387, 783 
Curran, P. A., Evans, P. A., de Pasquale, M., Page, M. J., & van der 
Horst, A. J. 2010, ApJ, 716, L135 
Dai, Z. G. & Cheng, K. S. 2001, ApJ, 558, L109 
Dai, Z. G. & Lu, T. 1998, Physical Review Letters, 81, 4301 
Dai, Z. G. & Lu, T. 1999, ApJ, 519, L155 
Daigne, F. & Mochkovitch, R. 1998, MNRAS, 296, 275 
Dainotti, M. G., Cardone, V. F., & Capozziello, S. 2008, MNRAS, 
391, L79 

De Colle, F., Ramirez-Ruiz, E., Granot, J., & Lopez-Camara, D. 2012, 
ApJ, 751, 57 

Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2009, MNRAS, 
397, 1177 

Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2007, A&A, 469, 
379 

Frail, D. A., Kulkarni, S. R., Sari, R., et al. 2001, ApJ, 562, L55 
Frail, D. A., Waxman, E., & Kulkarni, S. R. 2000, ApJ, 537, 191 
Friis, M., De Cia, A., Kriihler, T., et al. 2015, MNRAS, 451, 167 
Gao, H., Lei, W.-H., Zou, Y.-C., Wu, X.-F., & Zhang, B. 2013, New 
Astronomy Reviews, 57, 141 

Granot, J., Konigl, A., & Piran, T. 2006, MNRAS, 370, 1946 
Granot, J. & Piran, T. 2012, MNRAS, 421, 570 
Granot, J., Ramirez-Ruiz, E., & Loeb, A. 2005, ApJ, 618, 413 
Granot, J. & Sari, R. 2002, ApJ, 568, 820 

Greiner, J., Bornemann, W., Clemens, C., et al. 2008, PASP, 120, 405 
Greiner, J., Kriihler, T., Klose, S., et al. 2011, A&A, 526, A30 
Holland, S. T. & Pagani, C. 2012, GRB Coordinates Network, 13901, 
1 


12 








































Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 
440, 775 

Kirk, J. G., Guthmann, A. W., Gallant, Y. A., & Achterberg, A. 2000, 
ApJ, 542, 235 

Klebesadel, R. W., Strong, I. B., & Olson, R. A. 1973, ApJ, 182, L85 
Knust, F., Schady, P., & Greiner, J. 2012, GRB Coordinates Network, 
13891, 1 

Kobayashi, S., Piran, T., & Sari, R. 1997, ApJ, 490, 92 

Kriihler, T., Kiipcii Yolda§, A., Greiner, J., et al. 2008, ApJ, 685, 376 

Kumar, P. 1999, ApJ, 523, L113 

Kumar, P. & Barniol Duran, R. 2009, MNRAS, 400, L75 
Laskar, T., Zauderer, A., & Berger, E. 2012, GRB Coordinates Net¬ 
work, 13903, 1 

Lazzati, D. & Perna, R. 2007, MNRAS, 375, L46 
Leventis, K., van 2012MNRAS.427.1329L, H. J., Meliani, Z., & Wijers, 
R. A. M. J. 2012, MNRAS, 427, 1329 
Margutti, R., Zaninoni, E., Bernardini, M. G., et al. 2013, MNRAS, 
428, 729 

Medvedev, M. V. 2006, ApJ, 651, L9 
Meszaros, P. & Rees, M. J. 1997, ApJ, 476, 232 
Morsony, B. J., Lazzati, D., & Begelman, M. C. 2007, ApJ, 665, 569 
Nousek, J. A., Kouveliotou, C., Grupe, D., et al. 2006, ApJ, 642, 389 
Pagani, C., Barthelmy, S. D., Baumgartner, W. H., et al. 2012, GRB 
Coordinates Network, 13886, 1 

Page, K. L., Maselli, A., Mangano, V., et al. 2012, GRB Coordinates 
Network, 13892, 1 

Panaitescu, A. 2005, MNRAS, 363, 1409 
Panaitescu, A. & Kumar, P. 2000, ApJ, 543, 66 
Panaitescu, A. & Kumar, P. 2002, ApJ, 571, 779 
Pandey, S. B., Sahu, D. K., Resmi, L., et al. 2003, Bulletin of the 
Astronomical Society of India, 31, 19 
Pei, Y. C. 1992, ApJ, 395, 130 

Racusin, J. L., Liang, E. W., Burrows, D. N., et al. 2009, ApJ, 698, 
43 

Resmi, L., Ishwara-Chandra, C. H., Castro-Tirado, A. J., et al. 2005, 
A&A, 440, 477 

Rhoads, J. E. 1999, ApJ, 525, 737 

Santana, R., Barniol Duran, R., & Kumar, P. 2014, ApJ, 785, 29 
Sari, R. & Esin, A. A. 2001, ApJ, 548, 787 
Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 
Schlafly, E. F. & Finkbeiner, D. P. 2011, ApJ, 737, 103 
Schuller, F. 2012, in Society of Photo-Optical Instrumentation Engi¬ 
neers (SPIE) Conference Series, Vol. 8452, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, 1 
Schuller, F., Nord, M., Vlahakis, C., et al. 2010, BoA - The Bolome¬ 
ter Data Analysis Software User and Reference Manual, 4th edn., 
Max-Planck-Institut fur Radioastronomie, Argelander-Institut fur 
Astronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 
Siringo, G., Kreysa, E., Kovacs, A., et al. 2009, A&A, 497, 945 
Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 
Tanvir, N. R., Fynbo, J. P. U., Melandri, A., et al. 2012, GRB Coor¬ 
dinates Network, 13890, 1 

Tody, D. 1993, in Astronomical Society of the Pacific Conference Se¬ 
ries, Vol. 52, Astronomical Data Analysis Software and Systems II, 
ed. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes, 173 
van Eerten, H. 2014, MNRAS, 442, 3495 

van Eerten, H., Zhang, W., & MacFadyen, A. 2010, ApJ, 722, 235 
van Eerten, H. J. Sz MacFadyen, A. I. 2012, ApJ, 751, 155 
van Eerten, H. J. & Wijers, R. A. M. J. 2009, MNRAS, 394, 2164 
van Paradijs, J., Groot, P. J., Galama, T., et al. 1997, Nature, 386, 
686 

Waxman, E. & Meszaros, P. 2003, ApJ, 584, 390 
Wiersema, K., Covino, S., Toma, K., et al. 2014, Nature, 509, 201 
Wijers, R. A. M. J., Rees, M. J., & Meszaros, P. 1997, MNRAS, 288, 
L51 

Yost, S. A., Harrison, F. A., Sari, R., & Frail, D. A. 2003, ApJ, 597, 
459 

Zauderer, B., Laskar, T., & Berger, E. 2012, GRB Coordinates Net¬ 
work, 13900, 1 

Zhang, B., Fan, Y. Z., Dyks, J., et al. 2006, ApJ, 642, 354 
Zhang, B. & Meszaros., P. 2004, International Journal of Modern 
Physics A, 19, 2385 


13 



